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ORBITAL  TRANSFER  IN  MINIMUM  TIME 
V.  E.  BLEICK,*  U.  S.  Naval  Postgraduate  School 

1 .  Introduction.  The  problem  of  orbital  transfer  discussed 
here  is  that  of  scheduling  the  direction  p  of  constant  momentum 
thrust  of  a  rocket,  which  loses  mass  at  a  constant  rate,  so  that 
it  transfers  to  an  earth  satellite  orbit,  with  known  elements  of 
time,  position  and  velocity,  in  a  minimum  time  T  after  launching 
of  the  rocket.  The  launching  conditions  are  assumed  to  be  fixed. 
This  situation  is  illustrated  in  Figure  1  for  the  case  of  a  cir- 
cular orbit.  The  sector  angle  B  at  which  the  rocket  enters  orbit 
will  be  called  the  rendezvous  angle.  To  aid  the  discussion  imagin- 
ary physical  rendezvous  of  the  rocket  and  satellite  is  assumed  to 
occur  at  this  angle.  The  time  of  rocket  launch  to  achieve  actual 
physical  rendezvous  can  be  determined,  of  course,  only  after 
both  of  the  unknowns  T  and  B  have  been  found.  The  problem  is  set 
up  as  a  calculus  of  variations  problem  of  the  Lagrange  type,  and 
is  solved  by  an  iterative  process  in  which  an  initial  approxima- 
tion to  the  angle  B  is  estimated. 

A  non-rotating  Oxy  rectangular  coordinate  system  with  origin 
at  the  earth's  center  is  used.  The  coordinates  and  velocity  com- 
ponents of  the  rocket  and  target  satellite  are  denoted  by  x,y,u,v 
and  X,Y,U,V  respectively.  For  simplicity  the  equations  of  motion 
of  rocket  and  target  will  be  written  in  a  "non-dimensional"  form 
by  the  use  of  suitable  units.  The  unit  of  length  is  taken  as  the 
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Fig.  1 .  Orbital  transfer 


earth's  equatorial  radius,  R   -  20,925,000  feet.  The  unit  taken 
for  time  t  is  the  time  required  by  a  hypothetical  earth  satel- 
lite, in  equatorial,  circular,  vacuum,  sea  level  orbit,  to 

traverse  a  sector  of  one  radian.  This  unit  of  time  is  YR~7g   = 

e  ° 

13.459  minutes,  where  g  =  32.086  ft. /sec.   is  the  acceleration 
of  gravity  at  the  equator.  The  unit  of  velocity  is  then  the  speed 
of  this  hypothetical  satellite.  These  units  of  length  and  time 
will  always  be  understood,  unless  other,  more  conventional  units 
are  specifically  mentioned. 

2.  Statement  of  the  problem.  The  equations  of  motion  of  the 
rocket,  in  terms  of  the  specified  units  of  length  and  time,  are 

(1)        •  cp.  =  u   -   g      -  a  cos  p  =  0 
cp2=v-g2-a  sin  p  =  0 
cp   =  x  -  u  =  0 
^4  =  y  -  v  =  0 


3  /  3    2     2     2 

*  »  g2  =  -y/r  »  r    =  x    +  y 

where  m  is  the  constant  fraction  of  initial  gross  rocket  mass 


where  g1  =  -x/r  ,  g2  =  -y/r  ,  r   =  x   +  y  ,  a  =  cm/g(1  -  mt) , 


lost  per  unit  of  time,  c  is  the  constant  speed  of  the  emitted 
rocket  gases,  and  g  is  the  acceleration  of  gravity  at  the  equator. 
The  fixed  initial  conditions  of  the  rocket  trajectory  are  taken  as 

(2)        x(0)  =  0,   u(0)  =  V1  =  VQ  cos  6 

y(0)  =  1 ,   v(0)  =  V2  =  VQ  sin  6  . 
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The  terminal  point  of  the  rocket  trajectory  is  variable  with 

(3)     x(T)  =X(T),   y(T)  =  T(T) ,   u(T)  =  U(T) ,   v(T)  =  V(T) . 

It  is  also  assumed  that  the  rocket  thrust  is  turned  off  abruptly 
at  time  T.  This  discontinuity  will  lead  to  a  trivial  steering 
corner  in  the  calculus  of  variations  problem.  Note  that  (1)  and 

(3)  imply  that  U(T)  =  -[X/R3]T  =  g.,  (T)  and  V(T)  =  -[l/R3]T  =  g2(T) 
where  R  =  |T  X  +  J  T| . 

The  problem  is  to  choose  the  control  variable  p  to  effect 

/■i 

orbital  transfer  with   /~  dt   minimized,  and  to  determine  the 
corresponding  rocket  trajectory.  This  problem  is  equivalent  to 
the  Lagrange  calculus  of  variations  problem  of  requiring  the 
integral 

(4)  I   =  J0    (1    +   Xcp1    +  MCp2   +  ncp3    +   pcp4)dt 

to  be  stationary.  In  (4)  the  equations  (1)  are  regarded  as  con- 
straints with  X(t),  n(t)  ,  TT(t)  ,  p(t)  introduced  as  continua  of 
Lagrangian  multipliers  [1 ] .  If  the  time  coordinate  of  the  varied 
terminal  point  is  taken  as  T  +  AT ,  the  vanishing  first  variation.. 
[2]  of  I  is 

f  T 

(5)  61  =  J0    [X(6u  -  glx6x  -  g-|y6y  +  a  sin  P6P)  +  n(6x  "  6u) 

+  n(6v  -  g2x6x  "  g2y6y  "  a  cos  p6p)  +  p^y   "  6v) ldt 
f  T+AT 
+  JT    [1  +  X6u  +  n6v  -  (X  cos  p  +  [i    sin  p)6a]dt  =  0 
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where  the  finite  variation  6a,  6u,  6v  terms  in  the  integral  from 
T  to  T  +  AT,  created  by  thrust  termination  at  time  T  on  the  un- 
varied trajectory  and  at  T  +  AT  on  the  varied  trajectory,  cancel, 
On  integrating  by  parts  one  obtains 

(6)  61  =  [  X6u  +  (i6v  +  tt6x  +  p6y]T  -  L  [(X  +  tt)6u  +  (yi   +  p)6v 

+  (tt  +  glxX  +  g2xn)6x  +  (p  +  glyX  +  g2yn)6y 
+  a(n  cos  p  -  X  sin  p)6p]dt  +  AT  =  0. 

The  variations  of  the  dependent  coordinates  at  the  variable 
terminal  point  must  be  taken  [.2]  as 

(7)  6u(T)  =  (U  -  u)TAT  =  -(a  cos  p)TAT,   6x(T)  =  (X  -  x)TAT  =  0 
6v(T)  =  (V  -  v)TAT  =  -(a  sin  p)TAT,   6y(T)  =  (T  -  y)TAT  =  0. 

Substitution  of  (7)  into  (6),  and  application  of  the  fundamental 
lemma  of  the  calculus  of  variations  and  the  theory  of  ordinary 
extrema,  gives  the  Euler  equations 

(8)  X  +  tt  =  0 

|i  +  p  =  0 

A  +  SlxX  +  «2xU  =  ° 
P  +  g1y*  +  g2y^  =  o 
tan  p  =  (i/X 

and  the  transversality  condition 

[a(X  cos  p  +  |i  sin  p)  ]T  =  1. 
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The  first  four  homogeneous  equations  of  the  Euler  equations  (8) 

constitute  the  adjoint  system  [3]  of  the  system  of  variation 

equations 

(10)  6cp.  =  6cp2  =  6cp3  =  6cp.  =  0 

which  are  the  coefficients  of  X. ,  \if    rr,  p  in  (5)  equated  to  zero. 
The  adjoint  system  has  a  matrix  of  coefficients  which  is  the  neg- 
ative transpose  of  that  of  (10).  The  last  of  the  Euler  equations 
(8)  requires  that  the  control  variable  p  be  adjusted  so  that 
a  =  a(i  cos  p  +  j  sin  p) ,  which  is  proportional  to  the  rocket 
thrust,  is  continually  parallel  to  the  adjoint  vector  A  =  i  \  +  j  \i 
The  transversality  condition  (9) ,  which  may  be  written 
(a»A)   =  1  >  0,  requires  that  a  and  A  have  the  same  sense.  Since 
it  is  only  the  ratio  of  H  to  X  which  determines  p,  it  is  a  trivial 
matter  to  scale  them  to  satisfy  the  magnitude  requirement  of  (9). 

A  solution  of  the  problem  obtained  from  (1)  and  (8)  guarantees 
a  stationary  time  of  transfer.  The  nature  of  the  problem  is  such 
that  this  stationary  time  is  a  minimum  time. 

3.  Numerical  solution.  There  is  a  constructive  aspect  of  a 
modification  of  equation  (6),  first  used  by  Bliss  [4]  in  his  work*" 
on  differential  corrections  in  ballistics,  and  applied  recently 
by  Faulkner  [5]  in  an  iterative  fashion  in  optimum  control  prob- 
lems. To  find  the  desired  modification  of  (6),  assume  that  a 
solution  of  the  systems  (1)  and  (8)  has  been  obtained,  which  does 
not  necessarily  satisfy  the  terminal  conditions  (3).  Using  this 
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solution  and  holding  T  constant,  consider  the  variation  of  the 
vanishing  integral 
f  T 

(11)  J0     (Xcpn      4-    HCP2     +    TTCP3     +     pCp4)dt     =    0 

with  the  terminal  constraints  (3)  removed,  so  that  the  terminal 
variations  6u(T)  ,  6v(T)  ,  6x(T)  and  6(T)  become  free.  Since  X,  |~i , 
n,  p  satisfy  the  adjoint  system,  and  since  there  is  now  no  steer- 
ing corner  due  to  thrust  termination,  one  obtains 

TT 

(12)  [X6u   +  |i6v  +  n6x  +   p6y]m   =    /Q   a(-X    sin   p   +  M    cos    p)6pdt 

=  J0      A.(da/dp)6pdt 

where  a  =  a(i  cos  p  +  j  sin  p) .  Equation  (12),  which  is  the 
desired  modification  of  (6),  is  called  the  fundamental  formula 
by  Bliss  [4],  but  is  also  known  under  the  generic  name  of  Green's 
formula  [3].  By  the  use  of  (12)  it  is  possible  to  generate  the 
control  parameters  of  a  varied  trajectory  which,  hopefully, 
comes  closer  to  satisfying  the  desired  terminal  conditions  (3). 
To  do  this,  assume  that  the  adjoint  system  has  been  solved  to 
obtain  a  fundamental  set  of  four  linearly  independent  solutions 
given  by  the  rows  of 

(13)  Transpose  of  B(t)  =  [\.(t)  [x±{t)    TT.(t)  p±(t)  ]   i  =  1,2,3,4 

— » 

where  B(0)  =  I  is  the  identity  matrix.  The  solution  A  =  i  \    +  j*  |j 

of  the  adjoint  system,  required  to  satisfy  the  last  of  the  Euler 
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equations  (8),  is  taken  as  the  linear  combination 


(14)        X  =  X   +  IX   +  mX   +  nX 


4 


|J  =  |i1  +  l|i2  +  nm3  +  nn4 
so  that  the  control  angle  p  is  determined  by 

(15)    tan  p  =  (u1  +  m2  +  mu3  +  n^)/^  +  1X2  +  mX3  +  nX4) 
and  its  variation  by 


(16)    6p  =  [(Xn0^X0)6l  +  (X|i7-nXj6m  +  (^-^  )6n]/(  X^  +  u2) 


When  (13)  and  (16)  are  substituted  into  (12)  there  results  the 

system  of  Green's  formulae 

(17)      [6u  6v  6x  6y]TB(T)  =  [0   61   6m  6n]A 

where  the  elements  of  the  matrix  A  are 


(18) 


lj  =J0   a(^i  "  ^Xi)(X^j  "  ^j)<lt/(X2  +  M2)3/2. 


The  coordinates  of  the  terminal  point  of  the  varied  trajectory 
at  time  T  +  AT  may  be  taken  as  [u+Au   v+Av   x+Ax   y+Ay]„  where 

(19)  [Au   Av   Ax   Ay]T  =  [6u   6v   6x   5y]T  +  [u  v  x  y]TAT. 

In  an  effort  to  make  this  new  terminal  point  come  closer  to  satis- 
fying the  terminal  conditions  (3),  one  may  take 

(20)  [Au  Av  Ax  Ay]T  =  [U-u  V-v  X-x  Y-y]T  +  [U   V   X   Y]TAT. 
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Substitute  (19)  and  (20)  into  (17)  to  obtain 

(21)  [U-u  V-v  X-x  Y-y]TAT  +  [0  6l  6m  6n]AB~1 (T)  =  [U-u  V-v  X-x  Y-y]T 

i 

as  the  system  of  equations  for  the  determination  of  AT,  6l,  6m,  6n 
on  the  varied  trajectory.  The  Faulkner  [5]  scheme  for  the  numerical 
solution  of  optimum  control  problems  may  now  be  stated:  Make  an 
initial  guess  for  the  values  of  T,  1,  m,  n;  carry  out  a  simultan- 
eous numerical  integration  of  the  systems  (1),  (8)  and  (18)  using 
the  control  variable  of  (15);  solve  the  system  (21)  for  AT  and 
the  changes  in  the  control  parameters;  iterate  until  convergence 
is  obtained.  This  program  may  be  carried  out  in  a  matter  of  seconds 
on  modern  digital  computers. 

To  give  an  example  of  this  control  optimization  in  the  present 
problem  a  circular  satellite  orbit  was  assumed  of  radius  R  = 

|i  X  +  j  T|  =  1.075699,  corresponding  to  an  altitude  of  300  statute 
miles  above  sea  level.  For  this  orbit  |T  U  +  ;f  V|  =  1/Vr"  and 

| i  U  +  j  V|  =  1 /R  .  The  assumed  rocket  launching  velocity  was 
VQ  =  |?  V1  +  J  V2|  =  0.585402,  launch  angle  6  =  0.928084  and  ren- 
dezvous sector  B  =  0.153840  as  in  Figure  1.  Also  assumed  were 
c  =  10000.9  ft. /sec.  and  m  =  0.00360583  sec.-1  The  odd  appearance 
of  these  figures  is  related  to  the  difficulty,  explained  in  the 
next  section,  of  obtaining  the  initial  "guesses"  T  =  0.289725, 
1  =  -0.223125,  m  =  -29.9875  and  n  =  19.0847  .  With  this  input  the 
process  converged  in  four  iterations  to  the  seven  significant 
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figure  results  T  =  0.2894592,  1  =  0.1840054,  m  -  -108.94383, 
n  =  67.95886  and  B  =  0.1536015  .  Figure  2  shows  the  resulting 
trajectory  and  rocket  thrust  directions  for  equal  time  intervals, 
excepting  the  interval  terminating  in  transfer,  which  is  reduced 
by  one  half. 

4.  The  initial  guesses.  The  domain  of  convergence  of  the 
iteration  scheme  of  the  last  section  appears  to  be  rather  limited 
in  the  present  control  problem,  requiring  initial  guesses  for 
T,  1,  m,  n  which  make  [U-u  V-v  X-x  Y-y]™  small  in  (21).  In 
the  numerical  example  given  here,  where  the  angle  B  is  small  and 
radius  R  nearly  unity,  the  equations  (1)  can  be  linearized  and  an 
exact  solution  of  the  linearized  transfer  problem  used  to  supply 
the  input  guesses  required  to  solve  the  non-linear  problem.  The 
gravitational  terms  g..  and  g~  of  (1)  were  replaced  by  the  linear 
terms  in  their  Baylor  expansions  at  x  =  0,  y  =  1 .  The  systems  (1) 
and  (8)  then  become  uncoupled.  Their  solutions,  using  (2),  are 
(22) 


t 

()  a.  cos  p  sin  (t-w)dw  +  V..  sin  t 


x  ,  J0  a 

ft 

u  =  Jq   a  cos  p  cos  (t-w)dw  +  V-  cos  t 

Y2y  =  JQ  a  sin  p  sinh  Y2(t-w)dw  +  V2  sinh  fit   -  (1 /YD  cosh  f2t   +  3/fl 

ft 
v  =  J0  a  sin  p  cosh  Y~2(t-w)dw  +  V2  cosh  f2t  -  (l/f2)sinh  f2t 

tan  p  =  [l  cosh  lT2t  -  (n/f2)sinh  Y2t]/(cos  t  -  m  sin  t) . 
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It  will  be  of  no  avail  to  attempt  to  solve  (22)  and  (3)  for 

T,  1,  m,  n  by  the  Newton-Raphson  method,  since  the  Newton-Raphson 

equations  are  in  fact  the  system  (21)  of  the  poorly  convergent 

iterative  routine  of  the  last  section.  A  substitute  procedure 

of  solving  (22)  and  (3)  for  V- ,  V~  ,  m  and  B  was  used.  Vhen  the 

transfer  problem  has  been  solved  for  a  sufficient  number  of  such 

sets  ( V1 ,  Vp ,  m,  B)  ,  a  basis  will  be  at  hand  for  obtaining  desired 

rocket  launching  conditions  by  interpolation  or  extrapolation. 

It  can  now  be  revealed  that  the  numerical  example  of  the  last 

section  really  had  its  genesis  in  the  assumptions  B  =  rr/16,  8  =  tt/4  , 

VQ  =  0.65,  m  =  0.0025  sec."   and  c  =  10000  ft. /sec.   It  was 

estimated  that  T  =  0.3,  1  =  1.25,  m  =  -1.0  and  n  =  9.5  would 

satisfy  (22)  and  (3).  Using  a  cluster  of  five  closely  spaced 

points  around  (T,l,m,n),  a  single  application  of  regula  falsi  [6] 

was  made  in  the  hope  of  improving  these  values  of  T,  1,  m,  n.  The 

particular  cluster  chosen  gave  the  output  values  T  =  0.289725, 

1  =  -0.223125,  m  =  -29.9875,  n  =  19.0847,  which  will  be  recognized 

as  the  initial  guesses  of  the  last  section.  The  output  residuals  were 

[X-x,  Y-y,  U-u,  V-v]T  =  [0.041389,  -0.020160,  0.212093,  -0.221912]. 

These  residuals  were  then  processed  by  the  linear  system 

/T 
(23)  (X-x)T  =  AV1  sin  T  -  T(T)dB  +  AcL  a  cos  p  sin  (T-w)dw/c 

(U-u)T  =  AV1  cos  T  -  V(T)dB  +  AcL  a  cos  p  cos  (T-w)dw/c 

Y^(Y-y)T  =  AV2sinh  Y2T  +  Y2X(T)dB  +  Ac/Q  a  sin  p  sinh  Y2(T-w)dw/c 

(V-v)T  =  AV2cosh  f2T  +  U(T)dB  +  Acf^a  sin  p  cosh  f2(T-w)dw/c, 
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derived  from  (22)  and  (3)  for  the  case  of  a  circular  orbit,  to 
obtain  c  ■+•  Ac,  launching  conditions  V.  +  AV-  and  V2  +  AV2 ,  and  ren- 
dezvous sector  B  +  dB  which  will  give  small  residuals.  The  value 
of  c  +  Ac  turned  out  to  be  far  from  the  desired  10000  ft. /sec. 
Instead  of  changing  the  value  of  c,  the  transformation 

(24)     10000  In  [1  -  (m  +  Am)x]  =  (c  +  Ac)  In  (1  -  mT ) , 

which  leaves  the  integral  I ~   adt  invariant,  was  used  to  change  m 
and  cause  c  +  Ac  to  approach  10000  ft. /sec.  on  iterating  (24), 
(22)  and  (23).  Seven  iterations  produced  the  input  values  of  V. , 
V_ ,  c,  m  and  B  used  in  the  last  section. 

5.  The  problem  of  two  fixed  end  points.  The  problem  in  which 
the  conditions  at  both  ends  of  the  rocket  trajectory  are  fixed, 
or  the  problem  in  which  rocket  and  target  satellite  achieve 
actual  physical  rendezvous  at  a  specified  angle  B,  can  also  be 
solved  by  the  methods  presented  here,  but  with  greater  converg- 
ence difficulties. 

Acknowledgement.  I  am  indebted  to  my  colleague,  Prof.  Frank 
D.  Faulkner,  for  suggesting  the  problem  of  this  paper  and  for 
discussions  of  difficulties  encountered. 
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CDC-1604  ASSEMBLY  LANGUAGE  COMPUTER  PROGRAM 
Set  content  of  K100+1  ■  SLJ  0  L3440  ZRO  0  OOOOO 


KO 


K10 


K20 


K30 


REM 
REM 
REM 
REM 
REM 
SIL 
STA 
STQ 
LDA 
FDV 
SLJ 
STA 
FMU 
FDV 
STA 
LDA 
SLJ 
STA 
FMU 
STA 
LDA 
SIU 
SLJ 
STA 
FMU 
STA 
LDA 
FMU 
STA 
LAC 
FMU 
STA 
LAC 
FMU 
STA 
LAC 
FMU 
STA 
LDA 
FDV 
SLJ 
FMU 
STA 
LDA 
FMU 
FDV 
STA 
LDA 
STA 
ENI 
SLJ 
FMU 
STA 
LDA 
SLJ 
FMU 
STA 
ENI 
ENI 
LDA 
STA 
UP 
SLJ 
STA 


1 

0 
0 
0 
0 

k 

0 
0 
0 
0 
0 
H 
0 
0 
0 
0 

1* 

4 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
h 
0 
0 
0 
0 
0 
0 
0 
0 
4 
h 
0 
0 
0 

1+ 

0 
0 
3 
0 
3 
3 
3 
k 
0 


NT 

THETATEM 

DELTHETA 

UNITY 

R 

SQROOT 

V 

V 

R 

VSQDR 

B 

TRIG+70 

COSB 

R 

CAPYFIN 

B 

K20  +  6 

TRIG 

SINB 

R 

CAPXFIN 

COSB 

V 

CAPUFIN 

SINB 

V 

CAPVFIN 

COSB 

VSQDR 

CAPVDFN 

SINB 

VSQDR 

CAPUDFN 

REARTH 

GACCEL 

SQROOT 

MOOT 

OMEGA 

C 

MOOT 

GACCEL 

A 

THETATEM 

THETA 

0 

TRIG+70 

VSTART 

VI 

THETA 

TRIG 

VSTART 

V2 

23 

0 

EL1INIT 

ELI 

/-l 

70007 

EL1INIT 


COMPUTE  LINEARIZED  TRAJECTORY 

(Bl)  =  NUMBER  OF  TRAJECTORIES 

(BU)  =  NUMBER  OF  FALSI  ITERATIONS  LESS  ONE 

(A)  =  FIRST  DEPARTURE  ANGLE 

(Q)  =  OELTA  DEPARTURE  ANGLE 

STORE  NO.  OF  TRAJS. 

STORE  1ST  DEPART  ANGLE 

STORE  DELTA  DEPART  ANGLE 

RECIP.  ORBIT  RADIUS 
RET. JUMP  TO  SQUARE  ROOT 
STORE  TARGET  SPEED 


STORE  TARGET  ACCEL. 

TARGET  SECTOR  AT  RENDEVOUS 

RET. JUMP  TO  TRIGONOMETRIC  SR. 


TARGET  COOR.  AT  RENDEVOUS 
STORE  BU 

TARGET  COOR.  AT  RENDEVOUS 

TARGET  VEL.COMP.  AT  RENDEVOUS 

TARGET  VEL.COMP.  AT  RENDEVOUS 

TARGET  ACC.COMP.  AT  RENDEVOUS 

TARGET  ACC.COMP.  AT  RENDEVOUS 
EARTH  RADIUS  AT  EQUATOR 
ACCEL. GRAVITY  AT  EQUATOR 
RECIPROCAL  SHULER  FREQ. 

NON-DIMEN.  MASS  LOSS  RATE 
ROCKET  NOZZLE  VELOCITY 
C  MDOT 

NON-DIMEN.  ACCELERATION 

BEGIN  OUTER  LOOP 

SET  DEPARTURE  ANGLE 

SET  NO.  OF  FALSI  ITERATES 

COSINE    THETA 

INITIAL  HORIZONTAL  VELOC. 

SINE  THETA 

INITIAL  VERTICAL  VELOC. 

SET  FALSI  COOR.  QUINTUPLE 


PRINT  COORS.  DECIMAL.   B3=0 
PARAMETER  WORD 
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K40 


K50 


K60 


K70 


KlOO 


ZRO 
ENI 
SLJ 
ZRO 
ZRO 
SLJ 
ZRO 
SLJ 
ZRO 
STA 
ZRO 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
STA 
INI 
INI 
ISK 
SLJ 
ENI 
SLJ 
SLS 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ENI 
SLJ 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
LDA 
STA 
LDA 
STA 
ISK 
SLJ 
ENI 
ENI 
SLJ 
ZRO 
LDA 
STA 
UP 
SLJ 
SAU 
ZRO 
EXF 
SLJ 
ZRO 
ZRO 
ENI 
SLJ 
STA 


0  T5INIT 
b    0 

4  700C7 
0  ELI  INIT 
0  T5INIT 
4  XYUV 
0  0 

4  70007 
0  0 

0  DELXFIN 
0  DELVFIN 
0  DELXFIN 

5  DELX1 

0  DELYFIN 
5  DELY1 
0  DELUFIN 
5  DELU1 
0  DELVFIN 
5  DELV1 
5  5 
3  3 

3  23 
0  K40 
0  0 

4  MATRIX 
0  K50+2 
0  0 

0  0 

0  5 

0  5 

0  DELX1-1 

0  0 

0  0 

0  0 

0  DELXl+174 

0  0 

4  MATRIX 
0  K60 

0  0 

0  0 

0  4 

0  5 

0  DELXl+174 

0  5 

0  ELI 

0  4 

0  EL1+24 

EL2 

ELI 

DELX2-1 

DELX1-1 

23 

/-2 

20 

3 

XYUV 
0  0 

5  DELXFIN 
5  DELX5 

S  /-l 
4  70007 
0  EL5 
0  T5 
0  70 
4  70007 
0  EL5 
0  T5 

3  0 

4  70007 

0  DELXFIN 


SET  B5=0  FOR  K40+3 
PRINT  COORDINATES  OCTAL 
PARAMETER  WORD 

SET  FALSI  FUNC.  QUINTUPLE 

PRINT  FUNCTIONS  DECIMAL 

PARAMETER  WORD 


BEGIN  INNER  LOOP  OR  SLJ  0  K101  OR  L3440 
INVERT  DELXYUV  MATRIX 
SINGULAR  MATRIX  ALARM  HALT 


CODE  TO  INVERT  =  5 
L  =  5  ROWS  IN  DELXYUV  MATRIX 
ADD.  OF  DELXYUV  MATRIX  =  DELXl-1 
N=0.  WORD  NOT  USED 

M=0 

ADD.  OF  INVERSE  =  DELXl+174 

INVERSE  TIMES  LMNT  MATRIX 
WORD  NOT  USED 


MATRI X  MULT.  CODE  =  4 
L  =  5  ROWS  IN  INVERSE 
ADD.  OF  INVERSE  =  DELXl+174 
N=5  COLUMNS  IN  INVERSE 
ADDRESS  OF  LMNT  MATRIX  =  ELI 
M=4  COLS.  IN  MATRIX  PRODUCT 
ADDRESS  OF  PROD.  =  ELl+24 
SHIFT  COORS.  AND  FUNCS.   B3=0 


FOR  USE  IN  K70+2 

COMPUTE  NEW  DELXYUV  FUNCS. 


STORE  NEW  FUNCS. 


PRINT  NEW  COORS. 
PARAMETER  WORD 

CLEAR  ARITHMETIC 
PRINT  NEW  COORS. 
PARAMETER  WORD 


DEC.  WITH  PANEL 


ERRORS 
OCTAL 


CLEAR  B3 
PRINT  NEW 
PARAMETER 


FUNCS. 
WORD 


DECIMAL 
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Kl  10 


K120 
XYUV 


K130 


K140 


ZRO 

UP 

LDO 

LDA 

SLJ 

STQ 

ENI 

LDA 

STA 

UP 

LAC 

FMU 

FDV 

FAO 

STA 

SLJ 

ZRO 

SAU 

ZRO 

SLJ 

ZRO 

STA 

ZRO 

SLJ 

ZRO 

ZRO 

ZRO 

LDA 

FAD 

STA 

RSO 

AJP 

SLS 

SLJ 

LDA 

STA 

LAC 

FMU 

FAD 

STA 

FDV 

STA 

FMU 

FMU 

FMU 

FMU 

FMU 

FDV 

STA 

LDA 

SLJ 

FMU 

SLJ 

SLJ 

ZRO 

STA 

SLJ 

STA 

ENA 

STA 

ENI 

ENA 

STA 

UP 

ENI 

LDA 

SLJ 

STA 

FMU 

FS8 


0 

4 

0 

0 

4 

0 

1 

1 

1 

1 

0 

0 

0 

0 

0 

4 

0 

0 

0 

4 

0 

0 

0 

4 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

0 

4 

4 

0 

0 

4 

0 

0 

0 

1 

0 

1 

1 

0 
0 
4 
0 
0 
0 


DELVFIN 

K50-H 

ELS 

UNITY 

POLAR+130 

OQ 

2 

EL5 

EL 

/-) 

V 

TAUFIN 

R 

R 

S 

70007 

0 

CAPXFIN 

LAM 

70007 

0 

NT 

S 

70007 

0 

NT 

DELTHETA 

THETA 

DELTHETA 

THETATEM 

NT 

K20+5 

BEGIN 

0 

T1 

TAUFIN 

OMEGA 

TAUFIN 

UNITY 

STORl 

OMEGA 

ST0R2 

STOR2 

STOR2 

STOR2 

STOR2 

ERRCOF 

TAUFIN 

ST0R3 

STOR1 

LOG+66 

STOR3 

SQROOT 

SQROOT 

0 

DELTAU 

EQUAL 

DELT1 

4 

RUNKT1 

14 

0 

TAU1 

/-I 

0 

TAUFIN 

TRIG 

SINT 

VI 

CAPXFIN 


END  INNER  LOOP   (  OR  SLJ  0  L3440  ) 

DEPARTURE  THRUST  ANGLE 
MOVE  EL5,  EM5  AND  EN5 

(-V) 

(-VT) 

(-VT/R) 

B-( VT/R) 

S  =  INITIAL  TARGET  SECTOR 

ALARM  REENTRY  FROM  LOG  S.R. 

DECIMAL  DUMP  WITH  PANEL 

PARAMETER  WORD 

DECIMAL  DUMP 

PARAMETER  WORD 

OCTAL  DUMP 

PARAMETER  WORD 


SET  NEW  THETA 

END  OUTER  LOOP.  (SEE  K440) 
HALT.  START  RENDEVOUS  PROG. 
ENTER  DELXYUV  SUBROUTINE 


(-W) 
(-WT) 

ISTOR1 )  = 

(ST0R2)  = 

SQUARE 

CUBE 

FOURTH 

FIFTH 

MULT.  BY 


1  -  WT 
(1-WTJ/W 

- 1 5  .  B-  1  2 


1-WT 

LOG(l-WT)  TO  BASE  E 


COMPUTED  TAU  INCREMENT 
EQUALIZE  DELTAU  AND  SET  B2 
SET  EQUALIZED  DELT1 

SET  GILL  FOR  4"  INTEGRALS 


CLEAR  RUNKT1  AREA 


SIN  T 

VI  SINT 
VI  SINT  - 
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K150 


K160 


K170 


K200 


STA 

LDA 

SLJ 

STA 

FMU 

FSB 

STA 

LDA 

FMU 

SLJ 

ZRO 

STA 

FMU 

STA 

LDA 

FMU 

SLJ 

STA 

FMU 

FSB 

FAD 

STA 

LAC 

FMU 

FAD 

STA 

LDA 

FMU 

STA 

LAC 

FMU 

FAD 

FSB 

STA 

LDA 

FMU 

FAD 

SLJ 

ZRO 

STA 

LDA 

FDV 

STA 

FMU 

STA 

LDA 

FMU 

STA 

LDA 

FMU 

FMU 

STA 

LDA 

FMU 

FMU 

STA 

SLJ 

ZRO 

ZRO 

UP 

SLJ 

ENI 

SLJ 

ENI 

SLJ 

LDA 

FAD 

STA 

LDA 

FAD 


0 
0 
4 
0 
0 
0 
0 
0 
0 

s 

0 
0 
0 
0 
0 
4 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
3 
3 
0 
4 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
3 
0 
0 
0 
3 
0 
0 
4 
0 
0 
2 
0 
0 
4 
0 
0 
0 
0 
0 
0 
0 


DELXFIN 

TAUFIN 

TRIG+70 

COST 

VI 

CAPUFIN 

DELUFIN 

TAUFIN 

R00T2 

EXP+107 

0 

C0SHTR2 

ROOTHALF 

STOR1 

TAUFIN 

R00T2 

EXP+71 

SINHTR2 

V2 

STOR1 

THREDR2 

STOR1 

CAPYFIN 

ROOT2 

STOR1 

DELYFIN 

C0SHTR2 

V2 

STORl 

SINHTR2 

ROOTHALF 

STOR1 

CAPVFIN 

DELVFIN 

ELI 

ELI 

UNITY 

SQROOT 

0 

HYP 

A 

HYP 

STOR1 

SINT 

SINDOT 

STOR1 

COST 

COSDOT 

STOR1 

ELI 

SINHTR2 

SINHDOT 

STOR1 

ELI 

COSHTR2 

COSHDOT 

ALPHA 

RUNKT1 

DERIV1 

/+1 

/+3 

0 

ALPHA* 1 

0 

/-2 

DELXFIN 

SIN 

DELXFIN 

DELUFIN 

COS 


COS  T 

VI  COST 

VI  COST  -  U 


T  TIMES  SQ.  ROOT  2 


STORE  HYP.  COSINE 
C0SHTR2/SQ.R00T  2 


STORE  HYP.  SINE 
V2  SINHTR2 

ADD  (3/SQ.ROOT  2) 

(-Y) 


(STOR1 >  =  V2  C0SHTR2 
(-SINHTR2) 


ELI  SQUARED 
1  +  ELI  SQ. 


(STOR 1)  =  A/HYP  =  A  COS  P 

A  COSP  SINT 
A  COSP 

A  COSP  COST 

A/HYP 

A  SINP 

A  SINP  SINHTR2 

A/HYP 

A  SINP 

A  SINP  C0SHTR2 

SET  UP  GILL  ROUTINE 

PARAMETER  WORD 

JUMP  IF  B2  NOT  =  ZERO 
B2=0.  END  OF  INTEGRATION 
B2  NOT  =  ZERO 
INTEGRATE  AGAIN 

JUMP  TO  TEST  B2 
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K210 


DERIVI 


K220 


K230 


K240 


K250 


STA 
LDA 
FAD 
FDV 
STA 
LDA 
FAD 
STA 
SLJ 
LAC 
FMU 
FAD 
STA 
LDA 
FDV 
STA 
LDA 
FSB 
STA 
ENI 
SLJ 
STA 
LDA 
SLJ 
ZRO 
STA 
LDA 
FMU 
SLJ 
STA 
LOA 
FMU 
SLJ 
STA 
LDA 
SLJ 
ZRO 
FMU 
STA 
LDA 
SLJ 
FSB 
STA 
FMU 
STA 
LDA 
FMU 
SLJ 
ZRO 
FMU 
FDV 
STA 
LDA 
FMU 
SLJ 
FMU 
FS8 
STA 
FMU 
FAD 
SLJ 
STA 
LDA 
FDV 
STA 
LDA 
FOV 
STA 
FMU 
FMU 


0  DELUFIN 
0  DELYFIN 
0  SINH 
0  R00T2 
0  DELYFIN 
0  DELVFIN 
0  COSH 
0  DELVFIN 
0  XYUV 
0  TAUl 
0  OMEGA 
0  UNITY 
0  UNMINWT 
0  A 

0  UNMINWT 
0  ADOMINWT 
0  TAUFIN 
0  TAU1 
0  STOR1 
0  0 

4  TRIG 
0  SINT 
0  STOR1 
4  TRIG+70 
0  0 

0  COST 
0  STOR1 
0  R00T2 
4  EXP+71 
0  SINHTR2 
0  STOR1 
0  R00T2 
4  EXP+107 
0  C0SHTR2 
0  TAU1 
4  TRIG 
0  0 

3  EMI 

0  ST0R2 
0  TAUl 

4  TRIG  +  70 
0  ST0R2 

0  LAM 
0  LAM 
0  ST0R3 
0  TAUl 
0  R00T2 
4  EXP+71 
0  0 

3  ENI 

0  R00T2 
0  ST0R2 
0  TAU1 
0  ROOT2 

4  EXP+107 

3  ELI 

0  ST0R2 
0  MU 
0  MU 
0  STOR3 

4  SQROOT 
0  HYP 

0  MU 

0  HYP 

0  ST0R3 

0  LAM 

Q  HYP 

0  ST0R4 

0  SINT 

0  ADOMINWT 


TO  EXIT  OF  SUBROUTINE 

(-WT) 

1-WT 

A/(  1-WT) 

(STOR1 )  *  TAUFIN  -  TAUl 

(SINT)  =  SINE  OF  (STOR1) 


SIN  TAUl  *  -LAM3 


(STOR2)  =  -  EMI  LAM3 

COS  TAUl  -  LAM1 
LAM1  +  EMI  LAM3 
LAM  =  LAM1  +  EMI  LAM3 

(STOR3)  =  LAM  SQUARED 


SINH  TAU1R2 


(  -  ENI  MU4  ) 
(STOR2)  =  -  ENI  MU4 


COSH  TAU1R2  =  MU2 

ELI  MU2 

ELI  MU2  +  ENI  MU4 


LAM  SQ.  +  MU  SQ< 


(STOR3)  =  SINP 


(STOR4)  *  COSP 
COSP  SINT 


a  MU 
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STA  0 

SINDOT 

LDA  0 

ST0R4 

COSP 

FMU  0 

COST 

COSP  COST 

FMU  Q 
STA  0 

ADOMINWT 

COSDOT 

LDA  0 

ST0R3 

SINP 

FMU  0 

SINHTR2 

SINP  SINHTR2 

FMU  0 

ADOMINWT 

STA  0 

SINHDOT 

LDA  0 

STOR3 

SINP 

FMU  0 

COSHTR2 

SINP  C0SHTR2 

FMU  0 

ADOMINWT 

STA  0 

COSHDOT 

K260 

SLJ  0 

ALPHA+2 

EXIT  FROM  DERIV1 

ZRO  0 

0 

EQUAL 

SLJ  0 

0 

S.R.  TO  EQUALIZE  DELTAU 

LDA  0 

TAUFIN 

FDV  0 

DELTAU 

SLJ  4 

FIXIT 

CONVERT  DELTAU  TO  FIXED  PT. 

STA  0 

INTSIGN 

STORE  INTEGRAL  PART 

STQ  0 

FRACSIGN 

STORE  FRACTIONAL  PART 

AJP  1 

/  +  5 

JUMP  IF  A  NOT  =  ZERO 

QJP  0 

/  +  3 

A=0.  JUMP  IF  Q  =  0 

LDA  0 

TAUFIN 

A=0.  Q  NOT  =  0 

STA  0 

DELTFLSG 

ENI  2 

1 

SLJ  0 

/+12 

ENA  0 

0 

A  =  Q  =  0 

STA  0 

DELTFLSG 

K270 

ENI  2 

0 

SLJ  0 

/+10 

AJP  2 

/+1 

A  NOT  0.  JUMP  IF  A  GREATER  THAN  0 

SCM  0 

MASK 

A  LESS  THAN  0.  ABSO  INTEGER  TO  A 

INA  0 

1 

ABSO.  INTEGER  +  1  TO  A 

SAU  0 

/+1 

ENI  2 

0 

ABSO.  INTEGER  +  1  TO  82 

SCA  1 

2057 

LRS  0 

57 

ENA  1 

0 

LLS  0 

44 

STA  0 

FLOIPLON 

LDA  0 

TAUFIN 

FDV  0 

FLOIPLON 

STA  0 

DELTFLSG 

ENI  0 

0 

K300 

SIL  2 

FIXIPLON 

USED  IN  FINAL  PROG. 

LDA  0 

DELTFLSG 

SLJ  0 

EQUAL 

TO  EXIT  OF  S.R. 

ZRO  0 

0 

FIXIT 

BSS 

44 

Bl=5   B6=65302   P=460 

R 

DEC 

1.075698925  NON-DIMEN.  ORBIT  RADIUS 

V 

BSS 

1 

NON-DIMEN.  TARGET  SPEED 

VSQDR 

BSS 

1 

NON-DIMEN.  TARGET  ACCEL. 

B 

OCT 

1775622077325042   TARGET  SECTOR  AT  RENDEVOUS 

COSB 

BSS 

1 

SINB 

BSS 

1 

REARTH 

DEC 

20925000. 

EARTH  EQUATOR  RAD.  IN  FT. 

GACCEL 

DEC 

32.086 

EQUATOR  GRAVITY  IN  FT/SEC/SEC 

RMS 

BSS 

1 

OMEGA 

BSS 

1 

NON-DIMEN.  MASS  LOSS  RATE 

THETATEM 

BSS 

1 

VI 

BSS 

1 

NON-DIMEN.  INIT.  HORIZONTAL  VEL. 

V2 

BSS 

1 

NON-DIMEN.  INIT.  VERTICAL  VEL. 

ERRCOF 

DEC 

-15.B-12 

(-15/4096  DEC.)  USED  TO  FIND  DELTAU 

ROOTHALF 

OCT 

200055202363150C 

SINT 

BSS 

1 

COST 

BSS 

1 

SINHTR2 

BSS 

1 

COSHTR2 

BSS 

1 

UNMINWT 

BSS 

1 

STORE  FOR  (1-WT) 

PI/16 
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ELI  INIT 

EM  UNIT 

EN1  INIT 

Tl INIT 

EL2INIT 

EM2INIT 

EN2INIT 

T2INIT 

EL3INIT 

EM3INIT 

EN3INIT 

T3INIT 

EL4INIT 

EMU  IN  IT 

EN4INIT 

T4INIT 

EL5INIT 

EMS INIT 

EN5INIT 

T5INIT 

RUNKT1 

DELT1 

TAUl 

SINOOT 

SIN 

COSDOT 

COS 

SINHDOT 

SINH 

COSHDOT 

COSH 

K440 


EXP 

SQROOT 

TRIG 

AC.MINWT 

RC0T2 

TWO 

THREE 

MASK 

TENGRANO 

COMPRMS 

BEGIN 


LIO 


BSS 
OCT 
OCT 
OCT 
OCT 
OCT 
OCT 
OCT 
OCT 
OCT 
OCT 
OCT 
OCT 
OCT 
OCT 
OCT 
OCT 
OCT 
OCT 
OCT 
OCT 
BSS 
BSS 
BSS 
BSS 
BSS 
BSS 
BSS 
BSS 
BSS 
BSS 
BSS 
ENI 
ENI 
LDA 
STA 
UP 
LDA 
AJP 
SLS 
BSS 
BSS 
BSS 
BSS 
BSS 
OCT 
OCT 
OCT 
OCT 
DEC 
OCT 
REM 
ENI 
ENI 
ENI 
ENI 
LDA 
FDV 
STA 
ENI 
ENI 
ENI 
ENI 
ENA 
STA 
ENI 
ENI 
ENI 
ENA 
STA 


2 

200 

577 

200 

177 

200 

577 

200 

177 

200 

577 

200 

177 

200 

577 

200 

177 

200 

577 

200 

177 

1 

1 

1 

1 

2 

1 

2 

1 

2 

1 

3 

23 

0 

ELI 

ELI 

/-l 

NT 

K20 

BEG 

3 

1  17 

65 

1  16 

1 

200 

200 

200 

777 

100 

175 


150040 
637777 
446000 
6463  14 
150030 
637777 
446010 
646170 
150020 
637777 
446020 
646200 
150010 
637770 
446030 
646430 
150000 
637777 
446040 
646500 


0000000  DEC. VALUE  =  +1.251953125 

7777777  DEC. VALUE  =  -1.000000000 

0000000  DEC. VALUE  -  +9.500000000 

6314631  DEC. VALUE  =  +0.3000000000 

0000010  DEC. VALUE  *  +1.251464844 

7777777  DEC. VALUE  *  -1.000000000 

0000000  DEC. VALUE  ■  +9.503906250 

0000000  DEC. VALUE  ■  +0.2987060547 

0000000  DEC. VALUE  =  +1.250976562 

7777777  DEC. VALUE  =  -1.000000000 

0000000  DEC. VALUE  =  +9.507812500 

0000000  DEC. VALUE  =  +0.2988281250 

0000000  DEC. VALUE  =  +1.250488281 

7007777  DEC. VALUE  =  -1.000434756 

0000000  DEC. VALUE  =  +9.511718750 

0000000  DEC. VALUE  =  +0.3011474609 

0000000  DEC. VALUE  =  +1.250000000 

7077077  DEC. VALUE  =  -1.000006689 

0000000  DEC. VALUE  =  +9.515625000 

0000000  DEC. VALUE  =  +0.3017578125 


RESET  LMNT  INIT.  (SEE  K120) 


INIT 


+  5 

IN 


HALT.  START  RENDEVOUS  PROG. 


6  12 

0  0 
0  0 
0  0 
0  V 
0  R 
0  W 
0  0 
0  0 
0  0 
0  0 
0  35 
0  RUNKT 
0  0 

0  0 

1 8 

1  RUNKT+17 


B1  =  3 
Bl  =  6 
B1=7 
STORE 
155202363150 
240000000000 
260000000000 
777777777777 
CO. 

440000000000 
RENDE 
SET  N 


B6=65447  P=460 
86=65566  P=460 
86=65653  P=460 
FOR   A/U-WT) 


VOUS  PROG.  STARTS  HERE 
UMBER  OF  ITERATIONS 


V/R  =  W 

STORE  TARGET  ANGULAR  VEL, 


SET  GILL  S.R.  FOR  29  DIFF.EQUS. 


CLEAR  RUNKT+17  AREA 
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STA 

1 

RUNKT+20 

INI 

1 

2 

ISK 

1 

112 

SLJ 

0 

/~2 

LDA 

0 

UNITY 

STA 

0 

RUNKT+20 

SET 

LAM1(0)=1 

STA 

0 

RUNKT+37 

SET 

MU2(0)=1 

STA 

0 

RUNKT+56 

SET 

PI3(0)=1 

STA 

0 

RUNKT+75 

SET 

RH0U(0)=1 

STA 

0 

RUNKT+25 

SET 

PI1D0T(0)=1 

LAC 

0 

UNITY 

STA 

0 

RUNKT+47 

SET 

LAM3DOT(0)=-l 

STA 

0 

RUNKT+66 

SET 

MU400T(0)=-1 

LAC 

0 

TWO 

L20 

STA 

0 

RUNKT+44 

SET 

RHO2D0T(0)=-2 

ENI 

0 

0 

ENI 

0 

0 

ENI 

0 

0 

ENI 

0 

0 

ENA 

0 

0 

STA 

0 

RUNKT+2 

SET 

TAU  a  0 

STA 

0 

RUNKT+4 

SET 

X(0)=0 

LDA 

0 

VI 

STA 

0 

RUNKT+3 

SET 

XD0T(0)=V1 

STA 

0 

RUNKT+12 

SET 

U(0)=V1 

LDA 

0 

V2 

STA 

0 

RUNKT+6 

SET 

YDOT(0)=V2 

STA 

0 

RUNKT+15 

SET 

V(0)=V2 

LDA 

0 

UNITY 

STA 

0 

RUNKT+7 

SET 

Y(0)=1 

L30 

ENI 

0 

0 

SLJ 

0 

/+5 

BYPAS 

RECUR 

SIU 

6 

/+1 

SLS 

1 

/+1 

STOP 

ENI 

6 

0 

ENI 

1 

147 

LDA 

1 

SAVE 

STA 

1 

RUNKT 

UP 

1 

/-l 

ENI 

0 

0 

SLJ 

u 

EQUAL 

ZRO 

0 

0 

STA 

0 

RUNKT+1 

STORE 

SLJ 

4 

LAMMU 

RET, 

.JUMP  TO  OBTAIN 

SLJ 

0 

L70 

BYPAS 

ZRO 

0 

0 

LAMMU 

SLJ 

0 

0 

ENTRANCE     S.R. 

SIL 

1 

L60+3 

SAVE  B1 

SIU 

3 

L60+4 

SAVE  83 

LDA 

0 

RUNKT+20 

LOAD  LAM1 

STA 

0 

ST0R3 

LAM 

1  TO  ST0R3 

ENI 

0 

0 

ENI 

1 

0 

SET 

B1=0 

ENI 

3 

0 

SET 

B3=0 

LDA 

1 

EL 

FMU 

3 

RUNKT+34 

EL  (LA 

FAD 

0 

ST0R3 

STA 

0 

ST0R3 

INI 

3 

14 

ENI 

0 

0 

• 

ISK 

1 

2 

SLJ 

0 

/-3 

L50 

STA 

0 

LAM 

STORE 

FMU 

0 

LAM 

LAM 

SQUARED 

STA 

0 

STOR1 

LAM 

SQUARE  TO  STOR 

LDA 

0 

RUNKT+23 

LOAD 

STA 

0 

STOR3 

ENI 

3 

0 

LDA 

1 

EL 

FMU 

3 

RUNKT+37 

EL(MU2) 

- 

20  - 

160 


L70 


DERIV 


LI  10 


L120 


FAD 

STA 

INI 

ENI 

ISK 

SLJ 

STA 

FMU 

FAD 

STA 

SLJ 

ZRO 

STA 

FMU 

STA 

ENI 

ENI 

SLJ 

BSS 

LDA 

FMU 

FDV 

STA 

LDA 

FMU 

FDV 

FSB 

STA 

ENA 

SAL 

SLJ 

ENA 

SAL 

ENI 

SLJ 

LAC 

FMU 

FAD 

STA 

LDA 

FDV 

STA 

ENI 

LDA 

STA 

LDA 

STA 

LDA 

FMU 

STA 

LDA 

FMU 

FAD 

STA 

SLJ 

FMU 

STA 

FMU 

STA 

ENI 

SLJ 

LAC 

FDV 

STA 

LDA 

FMU 

FDV 

FAD 

STA 

INI 


0 

0 
3 
0 
1 

0 
0 
0 
0 
0 
4 
0 
0 
0 
0 

1 

3 

0 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
4 
0 
0 
0 
0 

1 

4 

1 

0 
0 
0 

1 

0 
0 

1 

1 


STOR3 

ST0R3 

14 

0 

2 

/-3 

MU 

MU 

STOR1 

HYP2 

SQROOT 

0 

HYP 

HYP2 

HYP3 

0 

0 

LAMMU 

3 

A 

LAM 

HYP 

RUNKT+1 1 

A 

MU 

HYP 

UNITY 

RUNKT+14 

/+2 

L220-4 

L  150-3 

ALPHA+2 

L22C-4 

0 

L220-3 

OMEGA 

RUNKT+2 

UNITY 

UNMINWT 

A 

UNMINWT 

ADOMINWT 

0 

RUNKT+12 

RUNKT+3 

RUNKT+15 

RUNKT+6 

RUNKT+4 

RUNKT+4 

ST0R2 

RUNKT+7 

RUNKT+7 

ST0R2 

STORl 

SOROOT 

STORl 

R32 

STORl 

R52 

0 

LAMMU 

RUNKT+4 

R32 

STORl 

ADOMINWT 

LAM 

HYP 

STORl 

RUNKT+1 1 

2 


STORE  MU 

MU  SQUARED 

LAM  SQ  +  MU  SQ 

HYP2  =  LAM  SQ  +  MU  SQ 


HYP  =  SQ  ROOT  OF  HYP2 

HYP3  =  HYP2  TIMES  HYP 

RESTORE  Bl 

RESTORE  B3 

TO  EXIT  OF  S-R. 

ENTER  FROM  L30+7 


SET  UDOTIO)  =  A(LAM)/HYP 


SET  VDOT(O) 

CURRENT  ADDRESS+2  TO  A 
SET  EXIT  FROM  DERIV.  PROG. 
JUMP  TO  SET  AIJDOT(O) 

REPAIR  EXIT  FROM  DERIV. PROG, 

BYPASS  DERIV.  PROG. 

l-WT) 

1-WT 

STORE  (1-WT) 

STORE   A/M-WT) 

SET  XDOT  =  U 

SET  YDOT  =  V 

LOAD  X 

X  SQUARED 

X  SQ  TO  ST0R2 

LOAD  Y 

Y  SQUARED 

X  SO  +  Y  SQ 

(X  SQ  +  Y  SQ)  TO  STORl 

STORE  RADIUS  CUBED 

STORE  RADIUS  FIFTH 

NOTE  Bl  SAVED  BY  LAMMU 
<-X)  TO  A 

(-X/R32)  TO  STORl 


SET  UDOT(T) 
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L130 


LI  40 


L150 


L160 


ENI 

0 

0 

ISK 

1 

5 

SLJ 

0 

L 12C-3 

ENI 

0 

0 

ENI 

3 

0 

LAC 

3 

RUNKT+26 

(-PI1 ) 

STA 

3 

RUNKT+17 

SET  LAM1DOT  =  -PI  1 

LAC 

3 

RUNKT+3  1 

(-RH0  1 ) 

STA 

3 

RUNKT+22 

SET  MU1DOT  ■  -RHOl 

LDA 

0 

RUNKT+7 

Y 

FMU 

0 

RUNKT+7 

Y  SQUARED 

FSB 

0 

STOR2 

Y  SQ  -  X  SO 

FSB 

0 

ST0R2 

Y  SO  -  21X  SQ) 

FMU 

3 

RUNKT+20 

(Y  SQ  -  2  X  SQ)LAMl 

STA 

0 

STOR1 

LAC 

0 

THREE 

MINUS  THREE 

FMU 

0 

RUNKT+4 

(-3X) 

FMU 

0 

RUNKT+7 

(-3XY) 

FMU 

3 

RUNKT+23 

(-3XY  MUD 

FAO 

0 

STOR1 

FDV 

0 

R52 

STA 

3 

RUNKT+25 

SET  PI  1DOT 

LAC 

0 

RUNKT+7 

<-Y) 

FMU 

0 

RUNKT+7 

{-  Y  SQUARED) 

FMU 

0 

TWO 

(-2  Y  SQUARED) 

FAD 

0 

STOR2 

X  SQ  -  2  Y  SQ 

FMU 

3 

RUNKT+23 

(X  SQ  -  2  Y  SQJMU1 

STA 

0 

STOR1 

LAC 

0 

THREE 

MINUS  THREE 

FMU 

0 

RUNKT+4 

(-3X) 

FMU 

0 

RUNKT+7 

(-3XY) 

FMU 

3 

RUNKT+20 

(-3XY  LAM1) 

FAO 

0 

STOR1 

FDV 

0 

R52 

STA 

3 

RUNKT+30 

SET  RHOIDOT 

INI 

3 

13 

ENI 

0 

0 

ISK 

3 

57 

SLJ 

0 

L13C-4 

LAC 

0 

OMEGA 

(-W).     ENTER  FROM  L7 

FMU 

0 

RUNKT+2 

(-WT) 

FAD 

0 

UNITY 

1-WT 

STA 

0 

UNMINWT 

LDA 

0 

A 

LOAD  INITIAL  ACCEL. 

FDV 

0 

UNMINWT 

A/( 1-WT) 

STA 

0 

ADOMINWT 

ENI 

0 

0 

LAC 

0 

MU 

FMU 

0 

RUNKT+20 

(-MU  LAM1) 

STA 

0 

STOR1 

LDA 

0 

LAM 

FMU 

0 

RUNKT+23 

LAM  MU1 

FAD 

0 

STOR1 

STA 

0 

Al 

Al  =  LAM  MU1  -  MU  LAMl 

LAC 

0 

MU 

FMU 

0 

RUNKT+34 

(-MU  LAM2) 

STA 

0 

STOR1 

LDA 

0 

LAM 

FMU 

0 

RUNKT+37 

LAM  MU2 

FAD 

0 

STOR1 

STA 

0 

A2 

A2  =  LAM  MU2  -  MU  LAM2 

FMU 

0 

Al 

Al  TIMES  A2 

FMU 

0 

ADOMINWT 

(A/1-WTJA1A2 

FDV 

0 

HYP3 

, 

STA 

0 

RUNKT+77 

SET  A12D0T 

LAC 

0 

MU 

FMU 

0 

RUNKT+50 

(-MU  LAM3) 

STA 

0 

STOR1 

LDA 

0 

LAM 

FMU 

0 

RUNKT+53 

LAM  MU3 
-  22  - 

FAD  0  ST0R1 

STA  0  A3 

FMU  0  Al 

FMU  0  ADOMINWT 

FDV  0  HYP3 

STA  0  RUNKT+102 

LAC  0  MU 

L170      FMU  0  RUNKT+6U 
STA  0  STOR1 
LDA  0  LAM 
FMU  0  RUNKT+67 
FAD  0  STOR1 
STA  0  AU 
FMU  0  Al 
FMU  0  ADOMINWT 
FDV  0  HYP3 
STA  0  RUNKT+105 
LDA  0  A2 
FMU  0  A2 
FMU  0  ADOMINWT 
FDV  0  HYP3 
STA  0  RUNKT+1 10 
LDA  0  A2 

L200      FMU  0  A3 

FMU  0  ADOMINWT 

FDV  0  HYP3 

STA  0  RUNKT+1 13 

LDA  0  A2 

FMU  0  AU 

FMU  0  ADOMINWT 

FDV  0  HYP3 

STA  0  RUNKT+1 16 

LDA  0  A3 

FMU  0  A3 

FMU  0  ADOMINWT 

FDV  0  HYP3 

STA  0  RUNKT+121 

LDA  0  A3 

FMU  0  A4 

L210      FMU  0  ADOMINWT 
FDV  0  HYP3 
STA  0  RUNKT+12U 
LDA  0  AU 
FMU  0  AU 
FMU  0  ADOMINWT 
FDV  0  HYP3 
STA  0  RUNKT+127 
EN  I  0  0 
SLJ  0  ALPHA+2 
ENI  0  0 
ENI  1  1U7 
LDA  1  RUNKT 
STA  I  SAVE 
UP  1  /-I 
SLJ  U  ALPHA 

L220      ZRO  0  RUNKT 
ZRO  0  DERIV 
UP  2  /+1 
SLJ  0  /+3 
ENI  0  0 
SLJ  U  ALPHA+1 
ENI  0  0 
SLJ  0  /-2 
LDA  0  CAPXFIN 
FSB  0  RUNKT+U 
STA  0  C3 
FMU  0  C3 
STA  0  ST0R2 
LDA  0  CAPYFIN 
FSB  0  RUNKT+7 


A3  =  LAM  MU3  -  MU  LAM3 
Al  TIMES  A3 


SET  A  13D0T 
(-MU  LAMU) 


LAM  MUU 

AU  =  LAM  MUU  -  MU  LAMU 
Al  TIMES  AU 


SET  A1UD0T 
A2  SQUARED 

SET  A22D0T 
A2  TIMES  A3 

SET  A23D0T 
A2  TIMES  AU 

SET  A2UD0T 
A3  SQUARED 

SET  A33D0T 
A3  TIMES  AU 

SET  A3U00T 
AU  SQUARED 

SET  AUUDOT 

END  OF  DERIV  PROG. 
ENTER  FROM  L70+7 

STORE  INIT.  CONDITIONS 
FOR  NEXT  ITERATE 

SET  UP  GILL  ROUTINE 
PARAMETER  WORD 

JUMP  IF  B2  NOT  =  ZERO 
B2=0.  JUMP  TO  RMS  COMP. 
B2  NOT  =  ZERO 
INTEGRATE  AGAIN 

JUMP  TO  TEST  B2 

BEGIN  R.M.S.  COMPUTATION 

CAPXFIN  -  XFIN 


CAPYFIN  -  YFIN 
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L230 


L240 


L250 


L260 


L270 


STA 
FMU 
FAD 
STA 
LDA 
FSB 
STA 
FMU 
FAD 
STA 
LDA 
FSB 
STA 
FMU 
FAD 
SLJ 
ZRO 
FDV 
STA 
ENI 
ENI 
LDA 
STA 
INI 
UP 
SLJ 
ZRO 
SLS 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
STA 
STA 
LDA 
STA 
STA 
LDA 
STA 
LDA 
STA 
STA 
LDA 
STA 
ENI 
ENI 
SLJ 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 


0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

1+ 

0 
0 
0 

1 

3 
3 

1 

3 
1 
4 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
3 

4 

0 
0 
0 
0 
0 
0 
0 
0 
0 


C4 

ST0R2 

ST0R2 

CAPUFIN 

RUNKT+12 

CI 

CI 

STOR2 

ST0R2 

CAPVFIN 

RUNKT+15 

C2 

C2 

STOR2 

SQROOT 

0 

TWO 

RMS 

17 

55 

RUNKT+20 

BB1  1 

-3 

/-l 

MATRIX 

0 

L250-3 

0 

0 

5 

4 

B811 

0 

0 

0 

BB1 1+120 

RUNKT+1 00 

BB1  1  +  100 

RUNKT+103 

BB1  1+101 

RUNKT+106 

BB1  1+102 

RUNKT+1  1  1 

BB1  1+103 

RUNKT+1 14 

BB1  1  +  104 

BB1 1+106 

RUNKT+1  17 

BB1  1+105 

BB1  1  +  11  1 

RUNKT+122 

BB1  1+107 

RUNKT+125 

BB1  1+1  10 

BB1 1+112 

RUNKT+1 30 

BB1  1+113 

0 

2 

MATRIX 

L270-2 

0 

0 

4 

4 

B81 1*120 

4 

8B1  1+100 

3 


CAPUFIN  -  UFIN 


CAPVFIN  -  VFIN 


STORE  ROOT  MEAN  SQUARE 

SET  4  BY  4  MATRIX 
LAMIMUIPIIRHOI  =  B 
(1=1,2,3,4)  AT  TAUFIN 

INVERT  THE  MATRIX 

SINGULAR  MATRIX  ALARM  HALT 


CODE  TO  INVERT  =  5 

L  =  4 

BB11  =  ADDRESS  OF  B  MATRIX 

N=0  NOT  USED  IN  INVERT 


M  =  0 
BB1 
A12 
SET 
A13 
SET 
A14 
SET 
A22 
SET 
A23 
SET 
SET 
A24 
SET 
SET 
A33 
SET 
A34 
SET 
SET 
A44 
SET 


NJOT  USED  IN  INVERT 
1+120  =  ADDRESS  OF  INVERSE 
TO  SET  4  BY  3  D  MATRIX 


01  1 
D12 
D13 
D21 

D22 

D31 

D23 
D41 


A12 

A13 

A14 

A22 

A23 
A32 

A24 
A42 


D32  =  A33 


D33 

D42 


A34 
A43 


D43  =  A44 


FOR  USE  6TH  WORD  HENCE 
MULT.  B  INVERSE  BY  D 
THIS  WORD  NOT  USED 


4  =  CODE  TO  MATRIC  MULTIPLY 

L  =  4  ROWS  IN  B  INVERSE 

BB1  1+120  =  ADDRESS  OF  B  INVERSE 

N  »  4  COLUMNS  IN  B  INVERSE 

BB1  1+100  =  ADDRESS  OF  D  MATRIX 

M  =  3  COLS.  IN  THE  MATRIX  PROD. 
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L300 


1310 


L320 


L330 


ZRO 
REM 
REM 
LDA 
STA 
UP 
ENI 
LDA 
STA 
UP 
ENI 
LDA 
STA 
UP 
LDA 
FSB 
STA 
LDA 
FSB 
STA 
LDA 
FSB 
STA 
LDA 
FSB 
STA 
SLJ 
STA 
ZRO 
ENI 
SLJ 
SLS 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ENI 
SLJ 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ENI 
SLJ 
STA 
ZRO 
ENI 
ENI 
ENI 
ENI 
ENI 
ENI 
ENI 
ENI 
UP 
SLJ 
SIU 
SLS 
ENI 


0  BB1 1+1 


BB1  1+12 

B42 

/-I 

2 

BB1  1  +  7 

B32 

/-l 

2 

BB1 1+4 

B22 

/-l 

RUNKT+1  1 

CAPUDFN 

BB1  1 

RUNKT+14 

CAPVDFN 

B21 

RUNKT+1 2 

CAPUFIN 

B31 

RUNKT+J5 

CAPVFIN 

B41 

L310 

BB1  1 

C4 

0 

MATRIX 

L310+1 

0 

0 

5 

4 

BB11 

0 

0 

0 

BB11+120 

0 

MATRIX 

L320-1 

0 

0 

4 

4 

BB1 1+120 

4 

CI 

1 

Cl  +  4 

0 

/+2 

BB11 

Cl  +  7 

0 

0 

0 

0 

0 

0 

0 

0 

/+1 

FINAL 

/+1 

/+! 

0 


BB1H-1  =  ADDRESS  OF  MATRIX  PROD. 
L273-306  SETS  UP  MATRIX  OF  SYSTEM  TO  BE 
SOLVED  FOR  DELTAUFIN,  VAREL ,  VAREM  AND  VAREN 


UDOTFIN 

UDOTFIN  -  CAPUDFN   (OR  PASS  IF  RENO.  POINT  FIXED 

VDOTFIN 

VDOTFIN  -  CAPVDFN   (OR  PASS  IF  REND.  POINT  FIXED 

j  Li 
UFIN 
UFIN  -  CAPUFIN   (OR  PASS  IF  REND.  POINT  F-IXED) 

VFIN 

VFIN  -  CAPVFIN   (OR  PASS  IF  REND.  POINT  FIXED) 

BYPASS  PRINT  OUT 
PARAMETER  WORD  TO  PRINT 


INVERT  SYSTEM  MATRIX 
SINGULAR  MATRIX  ALARM  HALT 


5  *  CODE  TO  INVERT 

L  =  4 

BBU  =  ADDRESS  OF  MATRIX 

N=0  NOT  USED 

M=0  NOT  USED 

BB1  1+120  =  ADDRESS  OF  INVERSE 

MULT.  INVERSE  BY  COL.  MATRIX< 
THIS  WORD  NOT  USED 


4  =  CODE  TO  MATRIC  MULTIPLY 

L  =  4  ROWS  IN  INVERSE 

BB1 1+120  =  ADDRESS  OF  INVERSE 

N  =  4  COLS.  IN  INVERSE 

CI  =  ADDRESS  OF  COLUMN  MATRIX 

M=l  COL.  IN  THE  COL.  MATRIX 

Cl+4=  ADD.  OF  COL.  MATRIX  RESULT 

BYPASS  PRINT  OUT 
PARAMETER  WORD  TO  PRINT 


JUMP  IF  B6  NOT  *  ZERO 
B6  ■  0.  JUMP  TO  FINAL  PROG- 
STOP  SWITCH 
RESET  B6 
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ENI 

1 

3 

LOA 

1 

CH4 

STA 

1 

VARTAU 

UP 

1 

/-l 

FAD 

0 

TAUFIN 

STA 

0 

SAVE+U3 

ENI 

0 

0 

L3U0 

FMU 

0 

W 

FAD 

0 

S 

STA 

0 

SPWT 

SLJ 

4 

TRIG+70 

STA 

0 

COSSPWT 

FMU 

0 

R 

STA 

0 

SAVE+136 

LDA 

0 

SPWT 

SLJ 

u 

TRIG 

ZRO 

0 

0 

STA 

0 

SINSPWT 

FMU 

0 

R 

STA 

0 

SAVE+135 

LDA 

0 

COSSPWT 

FMU 

0 

V 

STA 

0 

SAVE+137 

L350 

LAC 

0 

SINSPWT 

FMU 

0 

V 

STA 

0 

SAVE+11+O 

LAC 

0 

COSSPWT 

FMU 

0 

VSQDR 

STA 

0 

SAVE+142 

LAC 

0 

SINSPWT 

FMU 

0 

VSQDR 

STA 

0 

SAVE+141 

ENI 

0 

0 

ENI 

0 

0 

ENI 

1 

2 

LDA 

1 

Cl  +  5 

FAD 

1 

EL 

STA 

1 

EL 

UP 

1 

/-l 

L360 

SLJ 

4 

700C7 

ZRO 

0 

0 

SAU 

0 

MDOT 

ZRO 

0 

S 

EXF 

0 

70 

SLJ 

1+ 

70007 

STA 

0 

CAPXFIN 

ZRO 

0 

LAM 

SLJ 

0 

RECUR 

ZRO 

0 

0 

8SS 

3 

POLAR 

BSS 

200 

ELI 

8SS 

EMI 

8SS 

EN1 

BSS 

Tl 

BSS 

EL2 

BSS 

EM2 

BSS 

EN2 

BSS 

T2 

BSS 

EL3 

BSS 

EM3 

BSS 

EN3 

BSS 

T3 

BSS 

EL4 

BSS 

EMU 

BSS 

EN4 

BSS 

Jk 

BSS 

EL5 

BSS 

EM5 

BSS 

EN5 

BSS 

STORE  DELTAUFIN,VAREL,VAREM,VAREN 

NEW  TAUFIN  FOR  NEXT  ITERATE 

STORE  NEW  TAUFIN  (NOTE  L30+3 ) 

OR  SLJ  0  L350+5  IF  REND.  POINT  FIXED 


NEW  CAPYFIN 


NEW  CAPXFIN 


NEW  CAPUFIN 


NEW  CAPVFIN 


NEW  CAPVDFN 


NEW  CAPUDFN 


STORE  NEW  EL,  EM  AND  EN 

PRINT 

DECIMAL  DUMP  WITH  PANEL 

CLEAR  ARITH.  ERRORS 

PRINT 

DEC.  DUMP 

TO    NEXT    ITERATION 
Bl=33      B6=66370      P=U60 
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T5        BSS    1 

BSS    26         STORE  FOR  MATRIX  PRODUCT 

FINAL     ENI  1  147        COMPUTES  POINTS  ON  FINAL  TRAJ. 

LIL  2  FIXIPLON 

LOA  1  SAVE 

STA  1  RUNKT      RESET  INITIAL  CONDITIONS 

UP  1  /-I 

SLJ  t*  ALPHA      SET  UP  GILL  ROUTINE 

2R0  0  RUNKT      PARAMETER  WORD 

ZRO  0  DERIV 

ENI  0  0 

SLJ  4  /+4        RET. JUMP  TO  DUMP  SUBROUTINE 

UP  2  /+1        JUMP  IF  B2  NOT  =  ZERO 

SLJ  0  L70C+1     B2=0-  JUMP  TO  REMAINDER  OF  FINAL  PROG, 

L650      ENI  0  0 

SLJ  4  ALPHA+1    INTEGRATE  AGAIN 

ENI  Q  0 


SLJ  0  /-3 

SLJ  0  0  ENTRY  TO  DUMP  SUBROUTINE 

LDA  0  RUNKT+2    TAU 

ENI  0  0 

ENI  0  0 

STA  0  DUMP       (DUMP)  =  TAU 

LDA  0  RUNKT+4    X 

ENI  0  0 

ENI  0  0 

STA  0  DUMP+4     (DUMP+4)  =  X 

LDA  0  RUNKT+7    Y 

ENI  0  0 

ENI  0  0 
L660      STA  0  DUMP+5     (DUMP+5)  -  Y 

LDA  0  RUNKT+12   U 

ENI  0  0 

ENI  0  0 

STA  0  DUMP+6     (DUMP+6)  =  U 

LDA  0  RUNKT+15   V 

ENI  0  0 

ENI  0  0 

STA  0  DUMP+7     (DUMP+7)  =  V 

LDA  0  RUNKT+1  1   UOOT 

ENI  0  0 

ENI  0  0 

STA  0  DUMP+10    (DUMP+10)  =  UDOT 

LDA  0  RUNKT+14   VDOT 

ENI  0  0 

ENI  0  0 
L670      STA  0  DUMP+11    (DUMP+11)  =  VDOT 

LDA  0  RUNKT+12   U 

LDQ  0  RUNKT+15   V 

SLJ  4  POLAR+130 

ENI  0  0 

ENI  0  0 

STQ  0  DUMP+12    (DUMP+12)  =  ARCTAN(V/U) 

LDA  0  LAM 

LDQ  0  MU 

SLJ  4  POLAR+130 

ENI  0  0 

ENI  0  0 

STQ  0  DUMP+13    (DUMP+13)  =  ARCTAN(MU/LAM) 

SLJ  4  70007      PRINT 

STA  0  DUMP       PARAMETER  WORD 

ZRO  0  DUMP+13 
L700      SLJ  0  L650+2     TO  EXIT  OF  DUMP  S.R. 

ZRO  0  0 

LDA  0  C3         C3  =  CAPXFIN  -  XFIN 

ENI  0  0 

ENI  0  0 

STA  0  OUMP+14    (DUMP+14)  =  CAPXFIN  -  XFIN 

LDA  0  C4         C4  =  CAPYFIN  -  YFIN 


ENI  0  0 
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ENI 

0 

0 

STA 

0 

DUMP+15 

LDA 

0 

CI 

ENI 

0 

0 

ENI 

0 

0 

STA 

0 

DUMP+16 

LDA 

0 

C2 

ENI 

0 

0 

L710 

ENI 

0 

0 

STA 

0 

DUMP+17 

SLJ 

4 

700C7 

ZRO 

0 

0 

STA 

0 

DUMP+14 

ZRO 

0 

DUMP+17 

SLJ 

0 

L3570+5 

ZRO 

0 

0 

ST0R1 

BSS 

ST0R2 

BSS 

ST0R3 

BSS 

ST0R4 

BSS 

L720 

DEC 

•  • 

DELX1 

BSS 

DELY1 

BSS 

DELU1 

BSS 

DELV1 

BSS 

DEC 

i  • 

DELX2 

BSS 

DELY2 

BSS 

DELU2 

BSS 

DELV2 

BSS 

DEC 

i  • 

DELX3 

BSS 

DELY3 

BSS 

DELU3 

BSS 

DELV3 

BSS 

DEC 

•  • 

DELX4 

BSS 

DELY4 

BSS 

3ELU4 

BSS 

DELV4 

BSS 

DEC 

■  • 

DELX5 

BSS 

DELY5 

BSS 

DELU5 

BSS 

SAVE 

BSS 

150 

RUNKT 

BSS 

135 

CAPXFIN 

BSS 

CAPYFIN 

BSS 

CAPUFIN 

BSS 

CAPVFIN 

BSS 

CAPUDFN 

BSS 

CAPVDFN 

BSS 

TAUFIN 

BSS 

LAM 

BSS 

HYP 

BSS 

HYP3 

BSS 

MU 

BSS 

HYP2 

BSS 

UNITY 

DEC 

i  • 

W 

BSS 

OELTAU 

BSS 

SPWT 

BSS 

COSSPWT 

BSS 

SINSPWT 

BSS 

R32 

BSS 

R52 

BSS 

INTSIGN 

BSS 

FRACSIGM 

BSS 

FLOIPLON 

BSS 

FIXIPLON 

BSS 

(DUMP+15)  =  CAPYFIN  -  YF I N 
CI  =  CAPUFIN  -  UFIN 


(DUMP+16)  =  CAPUFIN  -  UFIN 
C2  =  CAPVFIN  -  VFIN 


(DUMP+17)  =  CAPVFIN  -  VFIN 
PRINT 

PARAMETER  WORD 

JUMP  TO  PRINT  B  FINAL 


ALSO  DELV5  STORE 
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THREDR2 

OCT 

2002417 

TREHALF 

DEC 

1.5 

DELTFLSG 

BSS 

NT 

BSS 

THETA 

BSS 

VSTART 

DEC 

.65 

DELTHETA 

RSS 

DELXFIN 

BSS 

DELYFIN 

BSS 

DELUFIN 

8SS 

DELVFIN 

BSS 

MOOT 

DEC 

.0025 

A 

BSS 

QQ 

BSS 

C 

DEC 

10000. 

EL 

BSS 

EM 

BSS 

EN 

BSS 

S 

BSS 

Al 

BSS 

A2 

BSS 

A3 

BSS 

A4 

BSS 

VARTAU 

BSS 

VAREL 

BSS 

VAREM 

BSS 

VAREN 

BSS 

ALPHA 

BSS 

134 

BSS 

14 

DUMP 

BSS 

20 

B81  1 

BSS 

B12 

BSS 

B13 

BSS 

B14 

BSS 

B21 

BSS 

B22 

BSS 

B23 

BSS 

B24 

BSS 

B31 

BSS 

B32 

BSS 

B3  3 

BSS 

B34 

BSS 

B41 

BSS 

B42 

BSS 

B4  3 

BSS 

B44 

BSS 

BSS 

120 

CI 

BSS 

C2 

BSS 

C3 

BSS 

C4 

BSS 

BSS 

4 

LOG 

BSS 

REM 
REM 

100 

L2000 

BSS 

REM 

1440 

L3440 

ENA 

0 

0 

ENI 

1 

17 

STA 

1 

BB11 

UP 

1 

/ 

LDA 

0 

TAUFIN 

SLJ 

4 

TRIG 

STA 

0 

BB1  1 

LDA 

0 

TAUFIN 

SLJ 

4 

TRIG+70 

ZRO 

0 

0 

STA 

0 

B21 

LDA 

0 

TAUFIN 

FMU 

0 

R00T2 

2002417416663160   THREE/R00T2 


NON-DIMEN.  INITIAL  SPEED 


FRACTIONAL  MASS  LOSS  PER  SEC, 


0000.     ROCKET  NOZZLE  VEL.  IN  FT. /SEC. 


INIT.  TARGET  SECTOR  ANGLE 


B1=30   B6=67340   P*460 


5  TIMES  L  SQUARED  CELLS 


STORE  FOR  MATRIX  PRODUCT 

B 1=12   B6=67700   P=460 

LOG  S.R.  ALARM  EXIT  TO  K110 

(67703)  =  750  65110  000  00000 

USNPGS    GEN.    DUMP         B1=32         P=2 1 

AID   TO    LINEARIZED    TRAJECTORY    PROG 


CLEAR    BB11    MATRIX    AREA 


(BBU  )    *    SINT 


(B21)    =    COST 
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L3450 


L3460 


L3470 


L3500 


L3510 


SLJ 

4 

EXP+71 

STA 

0 

B32 

LOA 

0 

TAUFIN 

FMU 

0 

R00T2 

SLJ 

4 

EXP+107 

STA 

0 

B42 

LDA 

0 

SIN 

FDV 

0 

C 

STA 

0 

813 

LDA 

0 

COS 

FDV 

0 

C 

STA 

0 

B23 

LDA 

0 

SINH 

FOV 

0 

C 

STA 

0 

B33 

LDA 

0 

COSH 

FDV 

0 

C 

STA 

0 

B43 

LAC 

0 

CAPYFIN 

STA 

0 

B14 

LAC 

0 

CAPVFIN 

STA 

0 

824 

SLJ 

0 

L357C+3 

STA 

0 

B34 

LDA 

0 

CAPUFIN 

STA 

0 

844 

LAC 

0 

DELXFIN 

STA 

0 

CI 

LAC 

0 

DELUFIN 

STA 

0 

C2 

LAC 

0 

DELYFIN 

FMU 

0 

R00T2 

STA 

0 

C3 

LAC 

0 

DELVFIN 

STA 

0 

C4 

ENA 

0 

/+2 

SAL 

0 

L320+4 

SLJ 

0 

L310 

ZRO 

0 

0 

ENA 

0 

L320+6 

SAL 

0 

L320+4 

LDA 

0 

Cl  +  4 

FAD 

0 

VI 

STA 

0 

STOR1 

FMU 

0 

STOR1 

STA 

0 

STOR2 

LDA 

0 

C2  +  4 

FAD 

0 

V2 

STA 

0 

ST0R3 

FMU 

0 

ST0R3 

FAD 

0 

ST0R2 

SLJ 

4 

SOROOT 

ZRO 

0 

0 

STA 

0 

VSTART1 

LDA 

0 

STOR1 

FDV 

0 

VSTART1 

SLJ 

4 

POLAR 

STA 

0 

THETA1 

LDA 

0 

C3+4 

FAD 

0 

C 

STA 

0 

EJECT 

LDA 

0 

C4  +  4 

FAD 

0 

B 

STA 

0 

SECTOR 

SLJ 

4 

70007 

STA 

0 

VSTART1 

ZRO 

0 

SECTOR 

SLJ 

4 

70007 

ZRO 

0 

0 

ZRO 

0 

VST ART  1 

(B32)  =  SINHTR2 


(B42)  =  C0SHTR2 

SINE  IMPULSE  INTEGRAL 


COSINE    IMPULSE    INTEGRAL 


SINH  IMPULSE  INTEGRAL 


COSH  IMPULSE  INTEGRAL 


JLNiyP    TC    PATCH 


CI  ■  ADD.  OF  MATRIX  COL 


CREATE  LINEAR  EQU.  ROUTINE  EXIT 
JUMP  TO  LINEAR  EQU.  ROUTINE 


REBUILD  RENDEVOUS  PROGRAM 
DEL  VI 

(ST0R1 )  -  NEW  VI 

(ST0R2)  =  NEW  VI  SQUARED 
DEL  V2 

(ST0R3)  =  NEW  V2 

NEW  VI  SQ.  +  NEW  V2  SQ. 

NEW  VSTART 

NEW  VI 

COSINE  NEW  THETA 

TO  INVERSE  COSINE  S.R. 

NEW  THETA 

DEL  C 

NEW  C 

(EJECT)  ■  NEW  C 

DIFFERENTIAL  OF  B 

NEW  B 

( SECTOR)  =  NEW  B 

PRINT  RESULTS  DECIMAL 

PARAMETER  WORD 

PRINT  RESULTS  OCTAL 

PARAMETER  WORD 
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L3520 


L3530 


13540 


VSTART1 
THETA1 
EJECT 
SECTOR 


ZRO 

ENI 

ENI 

LOA 

STA 

LDA 

STA 

LDA 

STA 

LDA 

STA 

SLJ 

ZRO 

STA 

FMU 

STA 

LDA 

SLJ 

ZRO 

STA 

FMU 

STA 

LDA 

FMU 

STA 

LAC 

FMU 

STA 

LAC 

FMU 

STA 

LAC 

FMU 

STA 

LDA 

SLJ 

ZRO 

FMU 

STA 

LDA 

SLJ 

FMU 

STA 

LDA 

FSB 

AJP 

SCM 

THS 

SLJ 

LDA 

FMU 

FDV 

STA 

LDQ 

SLJ 

BSS 

BSS 

BSS 

BSS 

SLS 

ZRO 

LAC 

FMU 

FAD 

AJP 

SLJ 

ZRO 

FMU 

FDV 

SLJ 


0  SECTOR 


0 
0 
0 
0 
0 
0 
0 
0 
0 
4 
0 
0 
0 
0 
0 
4 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
4 
0 
0 
0 
0 
4 
0 
0 
0 
0 
2 
0 
0 
0 
0 
0 
0 
0 
0 
0 


0 

0 

VSTART1 

VSTART 

THETA1 

THETA 

EJECT 

C 

SECTOR 

B 

TRIG+70 

0 

COSB 

R 

CAPYFIN 

B 

TRIG 

0 

SINB 

R 

CAPXFIN 

C0S8 

V 

CAPUFIN 

SINB 

V 

CAPVFIN 

COSB 

VSODR 

CAPVDFN 

SINB 

VSQDR 

CAPUDFN 

THETA 

TRIG+70 

0 

VSTART 

VI 

THETA 

TRIG 

VSTART 

V2 

C 

TENGRAND 

/+1 

MASK 

UNITY 

SECTOR+2 

C 

MOOT 

GACCEL 

A 

ELS 

K10C+2 

1 

1 

1 

1 

SECTOR+1 

0 

TAUFIN 

OMEGA 

UNITY 

SECTOR+1 

LOG+66 

0 

C 

TENGRAND 

EXP 


CONDITIONS 


RESUME  COMP.  OF  INIT, 
SET  NEW  VSTART 

SET  NEW  THETA 

SET  NEW  C 

SET  NEW  B 


TARGET  COOR.  AT  RENDEVOUS 


TARGET  COOR.  AT  RENDEVOUS 

TARGET  VEL.  COMP.  AT  RENDEVOUS 

TARGET  VEL.  COMP.  AT  RENDEVOUS 

TARGET  ACC.  COMP.  AT  RENDEVOUS 

TARGET  ACC.  COMP.  AT  RENDEVOUS 
COSINE  THETA 

INITIAL  HORIZONTAL  VEL. 

SINE  THETA 

INITIAL  VERTICAL  VEL. 

ROCKET  NOZZLE  VEL. 

NOZZLE  VEL.  DIFFERENCE 

JUMP  IF  DIFF.  IS  POS. 

ABSO.  VALUE  OF  DIFF. 

EXIT  IF  ONE  GREATER  THAN  DIFF. 

NO.  JUMP  TO  ADJUST  MOOT  AND  C 

YES.  LOAD  ROCKET  NOZZLE  VEL. 

NON-DIMEN.  ACCELERATION 

RETURN  TO  LIN.  TRAJ.  PROG. 
TEMP.  STORE  FOR  NEW  VSTART 
TEMP.  STORE  FOR  NEW  THETA 
TEMP.  STORE  FOR  NEW  C 
TEMP.  STORE  FOR  NEW  B 
HALT  IF  1-WT  IS  NEGATIVE 

(-T) 

(-WT)  . 

{  1-WT  ) 

HALT  IF  1-WT  IS  NEGATIVE 

LOG(l-WT)  TO  BASE  E 
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L3560 


L3570 


ZRO 
SCM 
FAD 
FOV 
STA 
LDA 
FDV 
SLJ 
ZRO 
FMU 
STA 
LDA 
STA 
ENI 
ENI 
ENI 
ENI 
ENI 
ENI 
ENI 
LDA 
FMU 
FDV 
STA 
SLJ 
LDA 
FMJ 
SLJ 
ZRO 
SLJ 
ZRO 
STA 
ZRO 
SLS 
ZRO 
END 


0 
0 
0 
0 
0 
0 
0 
4 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
4 
0 
0 
0 
0 
0 


0 

MASK 

UNITY 

TAUFIN 

OMEGA 

GACCEL 

REARTH 

SQROOT 

0 

OMEGA 

MOOT 

TENGRAND 

C 

0 

0 

0 

0 

0 

0 

0 

MOOT 

C 

GACCEL 

A 

K70 

CAPXFIN 

RC0T2 

L34d0+2 

0 

70007 

0 

UNITY 

SPWT 

L3570+7 

0 


NEW  OMEGA 

STORE  NEW  OMEGA 


SHULER  ANGULAR  FREQ. 

NEW  MOOT 
STORE  NEW  MDOT 
LOAD  NEW  C 
STORE  NEW  C 


STORE  NEW  A 

JUMP  TO  LIN.  TRAJ.  PROG. 

PATCH  FCR  L3460+1  L.I. 


ENTRY  FROM  L710+3 

PRINT  B  FINAL  DECIMAL 

PARAMETER  WORD 

FINAL  B  =  SPWT 

END  OF  NON-LINEAR  PROGRAM 
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TA7 
.U6 
no. 

34 

Bleick 

Orbital 
minimum  tii 

71 

IS  DEC  8f 
33  MAY  87 

58303 
transfer   ] 


1 


20b9t> 
->  Gi>  7  ^ 
3  1957 


58303 


Bleick 

Orbital  transfer  in 
minimum  time. 


genTA      7.U6  no. 34 

Orbital  transfer  in  minimum  time. 
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